Site-dependent NMR Spin-lattice Relaxation in the Superconducting State of an Iron 

Pnictide Superconductor 
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In a conventional superconductor, the spin-lattice relaxation rate on all nuclei should have the 
same temperature dependence below To- We performed 23 Na, 75 As, and 59 Co NMR studies on 
single crystals of NaFeo.95Coo.05 As, and found that spin-lattice relaxation rates show very different 
temperature dependent power-law behavior on three sites. We propose that such site-dependent 
behavior is due to the facts that the superconductor has two gaps of very different size. The power- 
law exponent of each nucleus is affected by the strength of the hyperfine coupling to the small gap. 
We also found that the large superconducting gap on the cobalt site is smaller than on other two 
sites. It suggests a local suppression of the superconducting gap on the dopant site. 

PACS numbers: 74.70.-b, 76.60.-k 



The Hebel-Slichter formulation is very successful in an- 
alyzing the NMR spin-lattice relaxation rate (SLRR) of 
superconductors^ 2 -. Below Tc, the scheme can be used 
to determine the symmetry and the amplitude of the su- 
perconducting gap. For example, the low-temperature 
SLRR usually shows a gaped behavior with a coherence 
peak in a s-wave superconductor, and a T 3 power law 
behavior in a d-wave superconductor. In the newly dis- 
covered high-Tp iron pnictide superconductors 3 - 6 , which 
are probably multiple-gap superconductors^—, however, 
the temperature dependence of the 75 As SLRR below T c 
varies dramatically with dopings. The SLRR of 75 As fol- 
lows a power-law temperature dependence by 1/Ti ~ T n , 
with the value of n varies from 6 to 1. 5 14-25 . It is unclear 
if the different value of n is caused by a change of the 
pairing symmetry with doping, a disorder scattering ef- 
fect in an 5± gap superconductor— ~—, or other unknown 
mechanism. 

So far, NMR studies in iron pnictides were mostly 
performed on 75 As. From the Hebel-Slichter formula- 
tion, the NMR SLRR of all nuclei on the same material 
should have the same temperature dependence, regard- 
less of s-wave or d-wave superconductivity. The underly- 
ing physics is that the SLRR is primarily determined by 
electron excitations, with the electron density of states 
(DOS) depending on the gap symmetry and the gap 
amplitude. In this letter, we present our NMR studies 
on high-quality NaFeo.95Coo.05As superconductors. The 
NMR is performed on three nuclei, 23 Na, 75 As and 59 Co. 
Below Tc, the SLRR show a power-law like behavior on 
all sites, but with very different power-law exponent n on 
three nuclei. Our analysis indicates that such non-scale 
behavior is well understood by two superconducting gaps 
but with different gap sizes. We further found that the 
large gap on the cobalt site is much smaller compared 
with other sites. Since cobalt serves as a dopant, such 
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FIG. 1: (color online) Left panel: The NMR spectrum (dotted 
lines) of 23 Na, 75 As, and 59 Co of a NaFeo.95Coo.05As single 
crystal at T =20 K (above T c ) and T =2 K (below Tc) with 
an applied field of 11.85 Tesla along the c-axis. The solid lines 
are the Lorentz function fitting. Right panel: A sketch of the 
NaFei-zCoz As unit cell. 



spatial variation of the gap value suggests a local sup- 
pression of gap amplitude on the dopant site. 

The growth of our electron-doped NaFeo.95Coo.05As 
superconducting single crystals has been reported 
previously^ 9 -. The same sample is used here, and the su- 
perconducting transition is indicated by the sharp tran- 
sition from the magnetization and resistivity at Tc ~ 18 
Our NMR measurements were performed on 23 Na, 
75 As, and 59 Co, with 6 Tesla and 11.85 Tesla magnetic 
field, and temperatures down to 1.5 K. 

As shown in Fig. [T] Na, Co and As atoms reside on 
three typical positions in the lattice 3 ^. Their typical 
NMR spectra are also shown in Fig.[TJ and all lineshapes 
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tional superconductorii^. Below T c , l/ 23 Ti, l/ 59 Ti 
and l/ 75 Ti show power-law-like temperature dependence 
(1/Ti w T n ), but the power law exponent are different 
with n ks 2 for 75 As, n w 1.7 for 59 Co, and n s» 1.0 for 
23 Na. The sequential decrease of n suggests more contri- 
butions from the low energy excitations. 

These SLRRs data cannot be fit by a single s-wave 
or a single d-wave function. Here we attempt to fit the 
data by a two-gap function, like in other iron arsenide 
superconductors 20 . The spin-lattice relaxation for a two- 
gap Fermi liquid superconductor is approximated by, 

poo 

1/T ln = 2j2 / (A} lf N*(E)) 2 f(E)(l-f(E))dE(2) 

i=l,2 J ° 

9 f°° 

Tm/Tu = — ]T / [d l s {E)ff{E){l - f(E)) dE (3) 



FIG. 2: (color online) The temperature dependent spin-lattice 
relaxation of 23 Na, 75 As and 59 Co of a NaFeo.95Coo.05As sin- 
gle crystal under a 11.85 Tesla magnetic field along the c-axis. 

are fit by Lorentz function. At T —20 K, a tempera- 
ture slightly above Tc, the full width of half maximum 
(FWHM) of the spectrum is about 8 kHz, 20 kHz and 34 
kHz for 23 Na, 75 As and 59 Co, respectively. The narrow 
linewidth at such a low temperature and a high field indi- 
cates that the quality of our sample is very high. Below 
Tc, the spectrum is significantly broadened due to the 
distribution of magnetic fields around the vortex core^i. 

Our spin-lattice relaxation is measured by an inversion 
pulse method. The SLRR in a solid follows 

wr^yj^,*!^), (1) 

where Ahf^q) is the hyperfine coupling constant, and 
Xj_ {q, ijj) is the imaginary part of the electronic dynami- 
cal susceptibility perpendicular to the magnetic field. 

The SLRR of a NaFeo.95Coo.05As single crystal un- 
der a 11.85 Tesla magnetic field is shown in Fig. [2] 
The 6 Tesla data (not shown in figure) gives identi- 
cal results, showing no field-dependent behavior. The 
SLRR results show that the hyperfine coupling orders 
as 59 Ahf > 75 Ahf > 23 Ahf. For comparison purpose, 
l/ 23 Ti and l/ 59 Ti are multiplied by a factor of 56 and 
0.11, respectively. Above Tc, the SLRRs of three sites 
show similar temperature dependence. From 80 K to 
200 K, \/T\T increases rapidly with temperature, and 
such behavior is known as a pseudogap-like phenomena 
reported in other iron pnictides^. Between 80 K and T c , 
the SLRR of all sites is weakly temperature dependent, 
which is an indication of spin fluctuations 3 ^. 

The Superconducting transition is shown by a 
drop of 1/Ti on all three sites. Similar to re- 
ports on other iron pnictides, no coherence peak is 
found in NaFeo.95Coo.05As, suggesting an unconven- 



where T\ n and T\ s are the spin-lattice relaxation time of 
the normal state and the superconducting state, respec- 
tively. Nq (i=l,2) is the normal state electron density 
of states (DOS) of each band, A l h * is the hyperfine cou- 
pling of the nucleus to each band, and f(E) is the Fermi- 
Dirac distribution function with f(E) = l/(e~ B/fcBT + l). 
Assuming a hyperbolic Fermi surface with Nq and Nq 
as constants, we define an effective normal sate DOS 
d = AifN&fDa ee AifNi/iAlfNl + A\ f Ni) with 
g?q + dp = 1, and an effective superconducting state 
DOS d l s {E) = A i hf N i s {E)/D with N l s as the DOS in 
each superconducting band. Then d l s is simplified as 
d l s = d Ng/NQ, with N^/Nq as a functions of two super- 
conducting gaps Ag and A|, depending on the detailed 
gap symmetry. In the end, Eq. [2] and Eq. [3] contain only 
four fitting parameters, D , d,Q, and A| . 

Since the system has spin fluctuations, we first fit the 
normal state SLRR data by 1/T^T ~ A/(T+Q) between 
Tc and 80 K. The 1/Ti data at all temperatures are then 
multiplied by (T + 0)/0 to remove the spin fluctuation 
effect so that the Fermi liquid assumption in Eq. [5]holds. 
Finally, we fit the data by Eq. [5] (above T c ) and Eq. [3] 
(below Tc) with the same parameters. 

TABLE I: The two-gap d-wave fitting parameters of the spin- 
lattice relaxation rate of 23 Na, 75 As, and 59 Co. 



site 


do 


Ag (meV) 


A| (meV) 


75 As 


0.78 


2.84±0.15 


0.26±0.05 


59 Co 


0.68 


1.81±0.15 


0.24±0.05 


23 Na 


0.57 


2.84±0.15 


0.20±0.05 



We first use a two-gap d-wave fitting, assuming the gap 
varies with angle by Aj(</>) = A l s cos(2(j)). The effective 
DOS, by averaging out the angle dependence, is given by 

4 = d lK(§-) for E > A* s , and 4 = d ^K(^) 

for E < Ac, where K(x) is an elliptic function. A 
good fitting is obtained for the three nuclei as shown in 
Fig.GDa), and the total effective DOS (= 4 + 4) below 
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FIG. 3: (color online) (a) The low-temperature 23 Na, 75 As, 
and 59 Co spin-lattice relaxation rate (symbols) and the fit- 
tings (solid lines) by the two-gap d-wave symmetry as de- 
scribed in the text, (b) The effective electron density of states 
on 23 Na, 75 As, and 59 Co sites obtained from the fittings. 



Tc are depicted in Fig. |3jb). The fitting parameters for 
each nucleus are listed in Table. Q] The close values of 
the large gap A<j and the small gap A| among all nuclei 
suggest that our fittings are physically reasonable. 

As shown in Fig. [3fb), the effective DOS is peaked at 
two energies corresponding to the gap value E = A l s 
and E — A|. The small gaps are similar for all sites 
with 2A|/fcsTc' ~ 0-3- The large gap Ag are similar on 
the As and the Na sites with 2A s /kgTc ~ 3.7, but is 
much smaller on the Co site with iA^/ksTc ~ 2.3. The 
smaller A^ on the cobalt site suggests a spatial suppres- 
sion of gap close to it. Since Co serves as a dopant in 
the lattice, the gap suppression is probably caused by a 
local scattering effect from doping. Similarly, a local gap 
suppression by dopant has been reported by STM studies 
in the cuprates 3 ^. 

In Table HI a systematic decrease of d^ (E = Ag) is 
seen in the order of 75 As, 59 Co, and 23 Na, which is re- 
versely demonstrated by the increase of the effective DOS 
close to small gap position (E = A|) in Fig. [3Jb). Phys- 
ically, it suggests that the 23 Na senses more low energy 
excitations than other sites, which explains the smaller 
power law exponent n of the SLRR on 23 Na. With d\ 
rewritten as d\ = 1/(1 + ( J 4 2 / /^ / )(iV ( 2 /iV 1 )), d\ is de- 
termined by the ratios of the actual DOS in two bands 
Nq/Nq and the ratio of the hyperfine coupling to two 
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FIG. 4: (color online) The increase of the NMR linewidth of 
23 Na and 75 As spectrum below Tc- Both linewidth are scaled 
to A~ 6 2 as shown by the vertical axis (see text), where X a b is 
the penetration depth. 



bands A\JA\p It is reasonable to assume that Nq/Nq 
does not vary with nucleus sites, and is determined by the 
two electron bands. A\JA^p however, could vary with 
nucleus site, since the hyperfine coupling is momentum 
dependent (see Eq. [1} , and the Fermi surfaces of different 
bands of the iron arsenides are separated in the Brillouin 
zone 3 ^. Then the final conclusion is reached that the 
smaller power-law exponent n of the SLRR is caused by 
a stronger local hyperfine coupling of the nucleus to the 
band with the small gap. 

Before going further, we confirm that our SLRR for all 
nuclei are intrinsic for several reasons. First, the spin re- 
covery for each nucleus has only one 7\ component below 
Tc- Second, the normal state \jT\ on all sites shows a 
similar temperature dependence, which is consistent with 
other electron-doped iron arsenide 3 ^. Finally and most 
importantly, we confirmed that 23 Na and 59 Co signals 
are excited from the same region around the vortex core 
within the penetration depth. The penetration depth is 
evaluated from the spectrum broadening below Tc due 
to the local distribution of fields around a vortex core 3 ^. 
For each nucleus, the penetration depth X a b is calculated 
by AB = Af/j n « O.O6O90 o /A 2 b below T c 31 , where <j> 
is the quantum flux, Af is the linwwidth, and j n is the 
gyromagnetic factor of the nucleus. As shown in Fig. [U 
the penetration depth approaches X a b ~0.3 /Ltm as T — s-0 
K from both 23 Na and 75 As measurements. Such pene- 
tration depth is typical for iron arsenide with a similar 
Tj£^ 

Therefore, our analysis suggests that the temperature 
dependence of the NMR SLRR below Tc is affected by 
two effects. First, the gap is locally suppressed by the 
disorder effect on the dopant site. Second, the different 
power-law exponent n of the SLRR on each nucleus is well 
understood in a superconductor with a large gap and a 
small gap. It is worthwhile to mention that ARPES may 
not be able to resolve the small gap. Although our fitting 
scheme is based on a d-wave symmetry, the conclusion is 
also valid on a two-gap superconductor with other sym- 
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metries. For example, we believe an anisotropic two-gap 
s-wave fitting should also work, if a similar electron DOS 
is obtained as that of the d-wave fitting. The proposed 
s± gap symmetry with strong scattering a 26 ' 27 ' 36 may also 
work. 

Our conclusion that the power-law behavior of the 
SLRR varies with local hyperfine couplings and local gap 
amplitude may shed light on 57 As NMR studies on other 
iron pnictides in two aspects. First, the local gap sup- 
pression from the dopant should be stronger with the in- 
crease of the doping. As a result, doping may cause more 
low energy excitations, and lower the power-law exponent 
as shown on the cobalt site. Second, with the increase of 
doping, the local hyperfine couplings Ahf{q) of 57 As on 
the Fermi surface may also vary with doping, because the 
Fermi surface changes dramatically with doping in iron 
arsenides 3 - 7 -. In both cases, the change of the power-law 
exponent n are not necessarily related to the change of 
the gap symmetry. 



In summary, our NMR study on NaFeo.95Coo.05As 
single crystals shows that the SLRR on 23 Na, 59 Co, 
and 75 As sites has very different temperature dependent 
power-law behavior below Tq- Our analysis suggests that 
such non-scale behavior is well described by the Hebel- 
Slichtcr formulation based on a two-gap superconductor 
with a large gap and a very small gap, regardless of the 
detailed gap symmetry. We also found that the large 
gap on the Co-dopant site is strongly suppressed, proba- 
bly due to the disorder effect. To our knowledge, this is 
the first report of a local gap suppression on the dopant 
site, and the effect may be verified by scanning probes 
such as STM. 
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